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Abstract 

We discuss the current use of the operator-product expansion in QCD calcula- 
tions. Treating the OPE as an expansion in inverse powers of an energy-squared 
variable (with possible exponential terms added), approximating the vacuum ex- 
pectation value of the operator product by several terms and assuming a bound on 
the remainder along the euclidean region, we observe how the bound varies with 
increasing deflection from the euclidean ray down to the cut (Minkowski region). 
We argue that the assumption that the remainder is constant for all angles in the 
cut complex plane down to the Minkowski region is not justified. 

Making specific assumptions on the properties of the expanded function, we ob- 
tain bounds on the remainder in explicit form and show that they are very sensitive 
both to the deflection angle and to the class of functions considered. The results 
obtained are discussed in connection with calculations of the coupling constant a s 
from the r decay. 

PACS numbers: 11.15.Tk, 12.38.Lg, 13.35.Dx 
July 1999 



1 e-mails: fischer@fzu.cz, Jan.Fischer@cern.ch 
2 e-mail: vrkoc@matsrv.math.cas.cz 



1 Introduction 



The operator-product expansion (OPE) [[T], ||, 

i f &x^ x A{x)B{<S) « C k {q)O k) (1) 
J k 

represents the product AB of two local operators as a combination, with c-number co- 
efficients, of the local operators Ok- Here, q is the total four-momentum of the system 
considered and q 2 = s = —Q 2 . The singularities of the product are contained in the 
coefficient functions Ck{q), which are ordered according to the increasing exponent k in 
s~ k . 

In local quantum field theory, the product Ai(x 1 )...A a (x a ) of two or several field oper- 
ators is singular for coinciding arguments, and the problem of defining it in a neighbour- 
hood is of fundamental importance. K. Wilson proposed that the operator product may 
be expanded in the form 

n 

A^x + x^.-.Aaix + x a ) = Y / fk(x\ ...,x a )O k (x) + Rn(x, x 1 ,..., x a ), (2) 

k=0 

where the remainder R n vanishes in the limit x J — > while the functions fk become 
singular or non-vanishing. He also generalized his hypothesis by assuming U that any 
operator product may be represented as a series 

oo 

A x (x + x 1 ).. .A a (x + x a ) = J2 fk(x\ ...,x a )O k (x), (3) 

k=0 

which is asymptotic in the sense that to every N there exists a k(N) such that the 
coefficients /^(pA 1 , ...,p\ a ) vanish faster than p N for all k > k(N). 

A rigorous formulation of the operator-product expansion within perturbation theory 
has been worked out by Zimmermann 0. The Wilson relation can be justified order 
by order in perturbation theory, and there is a well-defined algorithm allowing one to 
calculate the coefficient functions. Zimmermann, Wilson and Otterson H show how an 
operator product expansion can be derived from general principles, and find conditions 
under which the OPE gives complete information on the short-distance behaviour of 
operator products. 

Ferrara et al. [|J examine how the form of the operator-product expansion depends on 
the symmetry group of the theory. They find for instance that the covariance under the 
spinor group SU(2,2) can place significant restrictions on the structure of the expansion 
terms on the light cone. 

In momentum space, the Minkowski region lies along the cut Q 2 < carrying the 
spectrum of physical states, while in the euclidean region (Q 2 > 0) the expansion is, 
for large Q 2 , determined by short-distance dynamics, the separation of the large-distance 
contributions from the short-distance ones being well defined 0. Predictions in the 
Minkowski region are obtained by analytic continuation. 

The operator product expansion has been applied to various problems in quantum 
theory with varying degree of rigour. According to the problem considered, different 
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mathematical properties of the expansion have been proved or assumed (also the symbol 
~ in ([I]) is understood differently in different contexts). At a fixed perturbative order, 
one can express the operator product in terms of the Feynman diagrams of this order. 
There are reasons to believe that the large-momentum expansion has a non-vanishing 
convergence radius. Smirnov || (see also and references therein) proved asymptotic 
expansions of the renormalized Feynman amplitudes in the large- momentum (mass) limit, 
and found the corresponding operator expansions for the S-matrix and composite opera- 
tors. 

In quantum chromodynamics, a theory with a strong non-perturbative component, 
very little is known about the mathematical character of the operator-product expansion 
and its exact composition. In particular, it is not known whether terms exponential in 
the variable (Q 2 ) 1//2 /A, for instance terms of the form || 

exp(-C(g 2 ) 1 / 2 /A) (4) 

with C positive, have to be added to describe strong-interaction processes. There are 
reasons to believe that the series is divergent, but it is not known whether it is asymptotic 
to the function searched for and, if so, in what region of the Q 2 plane. We can expect 
that there is more chance to trust (|l|) at higher energies than at low energies; but the 
large-order behaviour of its terms is not known ||. 

In contrast to the lack of rigour, the operator-product expansion in QCD is very much 
needed as a tool for solving a number of practical problems, such as the semileptonic 
B-meson decay, heavy-light quark systems, heavy quarkonia and the Drell-Yan process. 
The inclusive decay hadronic widths are also expected to be calculable using the OPE 
and analyticity. 

To apply the operator product expansion in QCD, one is of course faced with the 
problem of extending it to non-perturbative dynamics. This issue has been discussed 
since the development of the QCD sum rules ||, where OPE is combined with analyticity 
and other non-perturbative aspects of QCD. The successful application of the QCD sum 
rule technique to many processes and effects is well known. 

Secondly, some information about the behaviour of the operator product expansion 
in the complex Q 2 plane away from euclidean region, along all rays passing through the 
origin, is necessary. The reason is that the Feynman graphs, through which the observables 
(including those describing high-energy effects) are expressed, contain integration over 
small momenta, where there is little chance that the OPE can be applied. Using the 
analyticity property, however, one can replace the low-energy integral by that along a 
circle of a sufficiently large radius in the complex plane (see Fig. 1). 

This, however, poses a new problem, that of finding conditions under which an expan- 
sion of the type 

/(1/Q 2 )^£a fc (g)/Q 2fc (5) 

k 

can be extended to the complex Q 2 plane, to be valid along all rays \Q 2 \ — > oo. Such 
an extension is a delicate problem requiring precise mathematical conditions, which are 
not known in the case of QCD. To make the problem well-defined, recourse to simplifying 
mathematical assumptions therefore seems necessary. 
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Figure 1: Integration contour in the s = —Q 2 plane, showing how analyticity is used to 
replace an integral along the cut by that along a circle. For the decay of the r lepton (see 
section 3), the radius is equal to its mass m T squared. 



While Wilson's operator product expansion is originally formulated in the Euclidean 
domain, its applications are mostly related to quantities of the Minkowski nature. Then, 
the assumption is usually adopted |L0| that the "convergence properties" of the operator- 
product expansion away from the euclidean ray Q 2 > are the same as those along it, 
except the points of the cut Q 2 < (the Minkowski region). Simultaneously, it is assumed 
that the truncation error (caused by approximating the operator product by several terms 
of the expansion) is independent of the direction in the complex Q 2 plane (see Fig. 1). 

We consider this assumption too a severe simplification, in addition technically moti- 
vated: indeed, it is hard to believe that the unknown discontinuity along the cut would 
have no influence on the behaviour of the function along rays that are near the cut. 

In the present paper we therefore look for a model scheme that would not a priori 
exclude the possibility that the bound on the truncation error (assumed originally in 
euclidean region) becomes looser with increasing deflection from the positive real semiaxis. 



The program of the paper was sketched in our previous papers, see [[Ll| and We 
discuss conditions under which an expansion of the type (^) can be extended to angles 
away from the euclidean semiaxis in the complex Q 2 plane, with the aim to find how a 
bound, originally assumed to be valid along the positive real semiaxis, develops when the 
deflection increases and the cut is approached. We propose (in section 2 of the present 
paper) a scheme that (i) has precise mathematical meaning, (ii) is free from the a priori 
assumption that the influence of the cut on the truncation error along a general ray 
can be neglected and, simultaneously, (Hi) tries to keep the model possibly close to real 
situations. To illustrate the physical interest of such a problem, we discuss in section 3 
the determination of the coupling constant from the r lepton hadronic width. Our results 
are presented and discussed in section 4. In the concluding section 5 we summarize our 
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results, and discuss several possibilities of refining the scheme to get closer to physical 
reality. The three Appendices contain the proof of the main theorem on which our results 
are based, and an example of a function that saturates the bound fl4"0|) obtained in the 
Appendix A. 

2 The operator-product expansion away from 
euclidean region 

Let F(s) be holomorphic in C, the complex s plane (s = — Q 2 ) cut along s G [0, oo] 
from which a bounded domain around the origin may be removed. Let the constants a&, 
k = 0, 1, 2, ...n — 1, and a positive number A n exist such that the following inequality 

n-1 

\F(s)-J2a k /(-s) k \<A n /\sr (6) 

k=0 

is satisfied for a positive integer n and all s < —b, with b being a positive number. (To 
avoid unnecessary complications, we take the lowest value of k to be zero, assuming that 
terms that are infinite at \s\ = oo have been removed.) The problem is what inequality 
(if any) will hold along rays in the complex s plane, away from the negative real semiaxis 
(euclidean region). 

As is natural to expect, the answer depends on additional assumptions imposed on 
the function F(s) = f(l/Q 2 ). When compared with the problem of the operator-product 
expansion in QCD, we make the following extra assumptions: 

1. The coefficient functions of local operators are perturbative series in powers 
of the QCD coupling constant a s (Q 2 ), which in turn is expanded in negative powers of 
ln(Q/A) with coefficient functions depending on lnln(Q/A), lnlnln(Q/A), etc. Here, A is 
Aqcd, the fundamental scale of quantum chromodynamics. It is sometimes instructive to 
consider the form of the power corrections in the case that the logarithmic Q 2 -dependence 
of the coefficient functions is neglected. We make this approximation for simplicity in the 
present paper, with the aim to examine the general case in a later publication. (In 
this approximation, the cut due to the logarithmic dependence of the Wilson coefficients 
disappears, but F(s) can still have a cut if the expansion has an infinite number of terms 
or if additional singular terms are added to the series, see below.) 

2. We assume that the inequality (|6]) is valid up to s = — oo , where the right-hand side 
of (P) vanishes at a high rate. This amounts to assuming that the n-th order remainder 

Rn(-1/S), 

n-1 

R n (z)=f(z)- y £a k z\ (7) 

k=0 

tends to zero for z — > as the n-th power of z for at least one value of n.f\ 

These two assumptions simplify our problem but they may move us farther from 
physics. We plan to refine the scheme in a subsequent paper; our principal intention 

3 We use the notation z — re lv — x + iy, z = — 1/s and f(z) = F(s). 
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here is to create a model which would be free from the conventional assumption that the 
truncation error is the same in all directions (except the cut). 

In phenomenological applications, the starting estimate of the truncation error in 
euclidean region is usually taken to be of the order of the first neglected term of the 
expansion. This term serves as a pragmatic guidance in estimating the size of the error; 
then, its value is conventionally extended to be the estimate of the truncation error along 
all rays passing through the origin. 

It is our aim to obtain a more realistic picture about how the error may develop when 
the cut is approached. We examine the high-energy properties of the expansion (|5|) along 
different rays in the complex plane. Assuming the bound @ for s < —b, we observe how 
it varies with increasing deflection of the ray from euclidean region. 

Note that we do not demand that the expansion (|5]) be convergent or asymptotic 
to the expanded function f{l/Q 2 ): our approach is more general and can be applied 
whenever the remainder R n [l/Q 2 ) obeys (|]) at least for one fixed value of n, i.e., tends 
to zero as 1/Q 2n or faster in the euclidean domain. This allows, under the conditions 
specified below, an analytic continuation of the truncation error from the euclidean to the 
Minkowski domain. If the starting inequality is known for several values of n, one can 
perform the continuation for each of them, term by term. 

Recent applications of the operator product expansion in QCD have focused on prob- 
lems dealing with quantities that are essentially related to the Minkowski domain, where 
the properties of the OPE are least known and may be completely different from those 
in the Euclidean domain. Inclusive decays of heavy flavours have been discussed. The 't 
Hooft model |13[ of two-dimensional QCD in the limit of many colours has been consid- 
ered, , with the aim to abstract general features that may survive in four- dimensional 
QCD. We refer the reader to the papers quoted for details. Here let us briefly discuss the 
example of the determination of a s {m^) from the r lepton hadronic width, to illustrate 
physical relevance of the problem. 



3 Determination of a s (m 2 T ) from the r lepton hadronic 
width 

The t lepton occupies a special position among all leptons, being the only lepton heavy 
enough to decay into hadrons. The r decay provides a unique chance to study hadronic 
weak interactions at moderate energies. There has been extensive interest in using mea- 
surements of its total hadronic decay width R T (normalized to the leptonic width), 

r[r~ — > v T + hadrons] 

to extract the renormalized strong-coupling parameter a s . This quantity possesses a num- 
ber of advantages compared with other QCD observables. It is expected to be calculable 
in QCD using analyticity and the operator product expansion. It is an inclusive quantity 
which has been calculated perturbatively to the order O(a^). The r mass, big as it is, is 
nevertheless below the threshold for charmed hadron production. 
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Starting from analyticity and the operator product expansion Braaten, Narison and 
Pich |TIJ used the measurements of the r decay rate to determine the QCD running 



coupling constant at the scale of the r mass m T . The ratio R T is represented in the form 

12tt / (l-s/m 2 T ) 2 (l + 2s/m 2 T )lmU(s) — , (9) 
Jo 

where IT(s) is a combination of correlators R^fy, corresponding to the two-point 

correlation functions for the vector V£ = ipj'j^tpi and axial vector = ^-7^75^ colour 
singlet massless quark currents, with coefficients given by the elements Kid and V ns of the 
Kobayashi-Maskawa matrix, the subscripts i, j =u,d,s denoting light quark flavours. For 
instance, 

U% v (s) = iJ d 4 xe^(0|T{^.(x),^(0)t}|0) = W-gTftUuAs). (10) 

The integral (§) cannot at present be calculated from QCD, because the hadronic functions 
are sensitive to the non-perturbative effects confining quarks in hadrons. But one can 
make use of the analyticity property of the correlating functions in the complex s-plane 
cut along the positive real semiaxis. This allows one to express (|9|) as a contour integral 
along the circle of radius m 2 : 

6m £ (l-s/m 2 T ) 2 (l + 2s/m 2 T )U(s)^ (11) 

(see Fig. 1, showing how analyticity is used to obtain R T from flll)). While in (|9|) the 



integration path < s < m 2 runs along the cut, the integration contour in flTT|) keeps 
distance from it (with the exception of one point, s = m 2 ., and its neighbourhood), thereby 
giving a justification for representing H(s) as the operator-product expansion over local 



operators (provided that the value of m T is large enough, see a discussion in Jl0|). 
In an analogous way, other weighted integrals (moments) of Imll(s), 

R 1 S/a(so) = J°° ds(l - s/s )\s/m 2 T r (12) 

have approximately been calculated within the framework of QCD. These moments can 
also be expressed as contour integrals analogous to (|iT|). 

It is to be expected that the error brought about by truncating the operator-product 
expansion of II(s) will be larger along rays that are closer to the cut s > 0. There is in 
particular a special danger that the integral (|ll|) receives essential contributions from an 
interval around s = m 2 , where the OPE has little chance appropriately to represent the 
function expanded. A fortunate circumstance is that the double zero of the kinematic 
factor (1 — s/m 2 ) 2 in the integrand suppresses the contribution from this dangerous seg- 
ment. But a quantitative analysis of this argument is, to the best of our knowledge, still 
lacking. Moreover, as is emphasized in experimental data are, because of the same 



factor (1 — s/m 2 ) I , very poor and statistically limited around this point. An explicit 
estimate of the error is therefore desirable. 
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It is our aim to develop a scheme allowing a quantitative analysis of this qualitative 
argument. The method of obtaining the value of ot s from R T was repeatedly criticised in 
the literature || |16| for various reasons, claiming that the theoretical uncertainty usually 
quoted is underestimated. Leaving aside uncertainties related to the truncation of the 
perturbation series, we focus on two aspects related to the truncation of the operator 
product expansion. These aspects are: 

1. The integral in flTT|) is carried out along the circle \s\ = m^. When H(s) in the 
integrand is approximated by the truncated OPE series, how does the remainder depend 
on the direction of the ray in the complex s plane? One rightly expects that the error will 
increase with s getting closer to the cut. But the quantitative aspect of this statement 
is unclear. How reliable is the expansion at those points \s\ — that are far from the 
euclidean interval, yet not quite close to the double zero at s = ? 

2. The terms exponential in (— s) 1 / 2 /A, possibly to be added to the operator product 
expansion, are small along the euclidean ray and invisible in the truncated series, but may 
become big when the ray in the s plane approaches the cut in the Minkowski domain. 
How does the effect of such additional terms depend on the deflection from euclidean 
region? 

In the subsequent sections, we discuss these problems by examining the mathematical 
model outlined in sec. 2. We find conditions for obtaining bounds on the truncation 
errors in the whole complex high-energy domain; it turns out that these bounds are larger 
(and die off slower) in the Minkowski domain than in the euclidean domain. 



4 Results and discussion 

According to the notation introduced in Sec. 2, f(z) represents a function holomorphic in 
a disk centred at the origin of the complex z plane cut along the negative real semiaxis. 
We consider two generic assumption schemes: 

Case I : Let us assume that f(z) has the form of the generalized Stieltjes integral (see 

/(*) = f r^-dt (is) 

Jo 1 + zt 

and that the moments _ 

a k = / t k p(t)dt (14) 
Jo 

exist for k = 0, 1, 2, n — 1, with n being a fixed positive integer and p(t) a real- valued 
function. Then, assuming a bound on the remainder (|7[) in euclidean region, we obtain a 
bound along all rays (see the formulae ([11]) and fliTf)). 

Case II : The integral representation (|T3"|) is not a dispersion relation, because the 
function p(t) has to vanish very rapidly for the moments a k (H) to exist. Beside this, 



non-power-like terms added to the operator-product expansion may violate the conditions 
(|j~3|), (0). We therefore build in subsection 4.2 a framework for more general situations, 
assuming an overall constant bound on R n (z). Then, using (||) on R n (z) along the eu- 
clidean ray, we obtain an angle-dependent bound along all rays (see the formula (|32|) 
or (PB])), and observe how it deteriorates when the ray deflects from euclidean to the 
Minkowski region. The result (|32"D, (|33"D is based on Theorem A of Appendix A. 
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In either case, the bounds obtained cannot be improved unless the respective class 
of functions is reduced to a smaller one. We show this by giving examples of functions 
saturating them in subsection 4.1 and in the Appendix C respectively. 



4.1 Case I : The remainder in the form of a Stieltjes integral 

Let us first assume that flTB| ) and flUf ) hold with p(t) nonnegative for t > 0. The remainders 



Rk(z) = (~z) k I " ^Tv/dt, (15) 
Jo 1 + zt 

with k = l,2,...,n — 1, are bounded, for z approaching zero along the positive real 
semiaxis, by 

\R k (z)\ = ~a k (z)\z\ k < a k \z\\ (16) 

where we use the notation 

a k (z) = r -f^j^t. (17) 
Jo \l + zt\ 

It is easily seen that the inequality (|16|) is valid also for z approaching zero along any ray 
lying in the right half of the complex plane because |i?fc(^)| < \A an d a k (z) < a k 
for Rez > 0. This, however, is not the case for Re z < 0, due to the presence of the cut 
in this halfplane. The following bound can be obtained for Rez < 0. We have, denoting 
z = re ltp , 

|1 + zt\ = (1 + 2rt cos p + r 2 t 2 ) 1/2 . (18) 

Considered as a function of t at t > and p fixed (vr/2 < \p\ < n), 1/|1 + zt\ has its 
maximum at t — — cosy?, where its value is 1/| sin^j. In this way we obtain (|1^) and 



\Rk(z)\ <a k \z\ k /\ sin <p\ (19) 

for Rez > and Rez < respectively. Comparing ( |T9"D with we see how the factor 
1/| sin (p\ makes the estimate looser when the ray gets closer to the cut, i.e., when p — > ±7r. 

The estimates become worse if the discontinuity along the cut is not positive definite. 
The corresponding bounds can be obtained by replacing p(t) with \p(t) \ in the derivation. 

Better bounds on |i?fc(z)| than ([IB]) and (|T^) cannot be obtained unless some special 
assumptions about pit) are made. To see this, choose pit) = po(t) such that f(z) has a 
pole, po(t) = c5(t — t ) with c > and t > 0, in which case f(z) = c/(l + zt ), a k = ctg, 
and \R k (z)\ = c|z| fc to7|l + ^o| fc - While the bound (|TBp for Rez > is saturated for t > 
tending to zero, ([19]) for Rez < is saturated for t = cos p. 

It is of interest to see how the bounds flTE| ) and (|H|) on the truncation error may affect 
the accuracy of the determination of the generic contour integral (0). As was discussed 
above, (Hd) and ( fT9] ) are not immediately applicable in QCD to estimate the truncation 



error, due to the simplifying assumptions of our model. But it is interesting to compare 
them with the conventional assumption that the error is constant in all directions of the 
complex z plane, which is a cruder approximation. To see this let us examine how the 
factor | siny9| _1 in (|l^) affects the estimate of the integrals of the form ([TT|) . 

As we have neglected the logarithmic dependence of the coefficient functions, the 



integral (|TT| ) can be evaluated trivially using Cauchy's residue theorem. Representing 
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the integration variable s in the form s = m^e lx , x — ft ~ ft an d replacing the function 
Il(s) by the bound ([K]), on Rk(z), we obtain the following bound on the integrated 
remainder: 



where 



and 



! Eu 



r -2k ( r2,l . r 2,l 

a k 67rm T (J Eu + I^ sin j . 

-3tt/2 



jl,m 
'Msin 



71-/2 



tt/2 

27T N 

3tt/2 / 



dx|l 



dx|l-e ix | / |l + 2e lx r/|sinx| 



2e lx 



(20) 



(21) 



[22] 



corresponds to the euclidean and the minkowskian halfplane respectively. 

The effect by which the discontinuity along the cut tells on the value of the truncation 
error can be seen when the sum 



rLm 



T l,m , j 
-'Eu + 1 



l,m 
Msin 



is compared with 
where 



n 



i,m 



jl,m 
'Eu 



Lm 



+ J Mi 



-'Mi 



tt/2 



+ 



2tt n 
3tt/2 



dx|l -e ix \ l \l + 2e lx | 



(23) 
(24) 
(25) 



in which the integrand does not contain the factor 1/| sin%|. The three integrals can be 
written in the form 

(26) 
(27) 



jl,m 
'Eu 



tt/2 



C"' m (x) d X , 



jl,m 
'Msin 



r /2 C l > m (x)/sm X d X , 
Jo 



and 



respectively, where 



jl,m 
i Mi 



tt/2 



C l > m ( x ) = 2'/ 2+m+1 (l-cosx) V2 (5/4 + cosx) m/2 . 



(2f 



(29) 



We see that (|20| ) is composed of two factors, Qna k m T (which depends on the fc-th 



moment and on the radius vn? T of the integration circle), and J 2 ' 1 = J 2 ' 



contains the factor 1/| sinx|). Certainly, 7 M ™ n > ^mTj which fact signals the zero of the 
denominator in (|15|). While /2' 1 is a divergent integral, turns out to be greater than 
I^ 1 by almost 24 per cent. This difference is 6.7 per cent for (l,m) = (2, 1) and further 
decreases with increasing I, but increases with increasing m at fixed I. Details can be seen 
in Table 1. 



En + J Msin (which 
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1 \ Tfl 


i 


o 
Z 


Q 
O 


A 


1 


1.237 


1.34 


1.45 


1.56 


2 


1.067 


1.104 


1.148 


1.194 


3 


1.026 


1.043 


1.064 


1.089 


4 


1.012 


1.020 


1.032 


1.046 



Table 1: The ratio I^™ / I{ m for different values of I and m. The ratio indicates how much 
the discontinuity along the cut affects the role of the truncation error in the integral I 1 ^" 1 
when compared with l[' m '. 



4.2 Case II : The remainder bounded by a constant 

For functions f(z) that do not satisfy the conditions of Case I, the following theorem may 
be useful. 

Let C be the plane of complex numbers. Let C(d, a), < a < ir, be the open segment 
of angle a, of the disk of radius d > centred at the origin. In other words, let C(d, a) 
represent z = re 1 *, \<p\ < a, < r < d (see Fig. 2). 

Theorem 1. Let a function f(z) be holomorphic in C(d,7c). Let the remainder R n (z), 
see (0), fulfill, for a fixed positive integer n, the following two inequalities: 

\R n {z)\<M (30) 

for complex z, z G C(d, it), and 

\R n {x)\<Ax n (31) 

for < x < d. Assume d = yM/A. Then, for every ip, \ip\ < tt, R n (z) satisfies the 
inequality 

\R n (re itp ) \ < M 1 ' 1 A 1 r in (32) 

for all z e C(d, 7r(l — 7)), where and 7 = 1 — \(p\/7r. 

This result is a special case of Theorem 2 (and Corollary 1), which is proved in the 
Appendix A. It gives an upper bound on R n (z) along every ray passing through the origin, 
the estimate becoming worse with increasing deflection from the positive real semiaxis, 
i.e., with the ray approaching the cut. 

Remark 1. The equality d = \j M/A is a special assumption saying that, for z ap- 
proaching the boundary point on the positive real semiaxis, z = d, the bounds (|30|) and 
(j31~l) coincide. We make this assumption for simplicity of the subsequent discussion; see 
Appendix A for the general case. 

Remark 2. The inequality (j32| ) can also be written in the form 

\R n (r^)\ <M(r/dy n , (33) 
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Figure 2: The region C(d, a) in the complex z plane. In this subsection, the case a = tt 
is considered. See the Appendix for the general case. 



which reveals how the bound deteriorates with increasing distance from the origin and/or 
increasing angle (i.e., when energy decreases and, respectively, when the cut is approached; 
see a discussion below). 

Remark 3. Compared with the conditions required in the previous subsection, no 
integral representations for f(z) or the coefficients are required in the Theorem. On 
the other hand, by (|3CiD a constant bound is imposed on f(z) in C(d, tt), which condition 
was not required in the Case I. In both cases, the resulting bound depends on the angle 
along which infinite energy is approached: In Case I, see ([16]) and (|I9D , the coefficient 
is angle-dependent and the exponent k is not, whereas in Case II, see (^) or (|3"3"D, the 
exponent of \z\ is angle-dependent, decreasing from n (along the positive real semiaxis) 
down to (along the cut). 

Remark 4- Since the resulting estimate (^) in C(d, tt) is considerably looser than (|3T| ) 
for < x < d (note that the exponent is changed from n in fl3l[) to in ( |32"D , < 7 < 1 , 
7 tending to zero near the cut), it is interesting to look for a function that saturates it. 
A set of functions saturating ([32] ) for different nonnegative integers n can be generated 
by using the function f(z) = exp{ia [ln(i/i: e~ 17r / 2 ) ] 2 } with a real. This example (see the 
Appendix C for details) shows that the bounds are optimal, in the sense that they cannot 
be improved within the class of functions considered. There might be physical reasons, 
on the other hand, to restrict oneself to a smaller class of functions, in which case an 
improvement of the bound would be possible. 

The resulting inequality (|32|) can be used to estimate, under the assumptions made, 
the error caused in (|TT|) due to approximating II (s) in (|ll|) by the first n terms of the 
expansion. We observe the following facts: 

1. To obtain a bound on the integral (|TT]) , we assume both (^) and (0). The 
inequality (RTF) alone (valid only in the euclidean region) is not sufficient for obtaining a 
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bound on the remainder R n (z) at complex z, unless ( ]30|) is simultaneously used. 

2. By Theorem 1, the bounds ( |3"0"D and ( J5T| ) are combined to create a third one, 
While (|30D is valid in the whole complex region C(<i, 7r) but is loose, (|3~l"D holds only on 
the segment < x < d but is far more restrictive on this interval. Theorem 1 combines 
them to produce (0), which holds in C(d, ir) and is restrictive, although less than (|3T]), 
becoming (|3TD and (|3~0"D on the positive and the negative real semiaxis respectively. The 
resulting bound M 1_7 A 7 r 7n on the right hand side of ( |32"D becomes ( |3T| ) and ( |30"D for 
7 = 1 (euclidean region) and 7 approaching zero (minkowskian region) respectively. 

3. For r = d, (K3) is no improvement of (Bp!). For r < d, fl3"2|) does imply an 



improvement of (|30"D thanks to (|3~lD , by means of the factor (r/d) 



„7n 



(A/ 'M) 7 on the 



right hand side. This factor is smaller than 1 for r < d and < ir, but approaches 1 
when \ip\ approaches 7r at fixed r, or when r approaches d at fixed tp. 

4. These results induce analogous relations between the corresponding integrals of the 
type (|), (|TT|) and ([12]). Let us introduce the notation p = r/d. Inserting Mp^ 1 '^^^, 
the right-hand side of (0), into (p7T|) instead of n(s), we find that the integral is bounded, 
due to (|30|), by 

Ml[' m (34) 

while (|3~2"D imposes the bound 
with 



MJ l ' m (n,p) 



J l > m (n,p)= / C"< m ( X ) p nx/ - d X . 
Jo 



(35) 
(36) 



on the same integral. 



n \ p 


0.05 


0.2 


0.4 


0.6 


0.8 


0.9 


1 


1 


0.18 


0.38 


0.57 


0.73 


0.87 


0.936 


1 


2 


0.051 


0.17 


0.34 


0.54 


0.76 


0.87 


1 


3 


0.020 


0.081 


0.21 


0.40 


0.66 


0.82 


1 


4 


0.009 


0.044 


0.13 


0.30 


0.58 


0.77 


1 



Table 2: Thanks to ([31]), (|30| ) improves into (|3^) and the corresponding integrated bound 
Ml l { m into MJ l ' m (n, p). The values of J 2,1 (n, p)//^' 1 shown in the table for different n 
and p illustrate the effect of Theorem 1 on the integrated bounds. Small values mean a 
big improvement, 1 means no improvement. 



As we have seen, ( |3lf) can be used to improve the original bound ([[(]) into (|32f) . This 
induces an improvement of the bounds on the corresponding integrals, changing Ml l { m 
into MJ l,m (n,p). The improvement is pronounced for small values of p (which means 
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that the disk C(d, n) is large), but becomes weak for p approaching 1, when the boundary 
circle of C(d,n) is approached. The values of the ratio J l,m (n, p)//{' m for (l,m) = (2,1) 
and some typical n and p are shown in Table 2. 

It is difficult to make a straightforward comparison of the two bounds, fl3"2] ) on one side 
and flTED, (^) on the other, because they have been derived under different conditions. 
Perhaps the most striking difference is that, in Case I, Rk{z) is holomorphic in the whole 
cut z plane, while in Case II it is holomorphic only in the cut disk C(d, it); in this way, 
singularities of any kind are allowed outside C(d, ir), arbitrarily near the boundary circle. 
The two assumption schemes therefore assign similar analyticity properties to f(z) if the 
disk C(d, it) is very large, i.e., for small values of p in fl36"l) . 



n \ p 


0.05 


0.2 


0.4 


0.6 


0.8 


0.9 


1 


1 


3.7 


1.9 


1.43 


1.21 


1.086 


1.040 


1 


2 


20 


4.1 


2.1 


1.49 


1.18 


1.081 


1 


3 


155 


10 


3.3 


1.9 


1.29 


1.13 


1 


4 


10 3 


27 


5.2 


2.3 


1.41 


1.17 


1 



Table 3: When extrapolated into the cut disk, (|31| ) deteriorates into (|32|) and the cor- 
responding integrated bound Mp n I l { m into MJ l ' m (n, p). The table shows the values of 
J 2,1 (n, p)/(p n /^' 1 ) for different n and p. Large values, which signal a big deterioration 
(due to a strong angle dependence of (|3"2|)), appear near the origin (where p is small) 
and/or at high expansion orders n. 



The constant M in (pOf ) affects the value of the integral bound ([35|) in combination with 
the damping factor the integrand of (|36|), which becomes small when n increases 

and/or p decreases. It is instructive to compare the resulting bound MJ l,m (n, p) with 
that based on the assumption that fl3"T|) preserves its form in the whole cut disk. A 



convenient way to measure this effect is to divide MJ l ' m (n, p) by Mp n /{' m , in which the 
angle-dependent factor 

p nxh i n 

the integrand is suppressed, the exponent nxj-n being 
replaced by its value along the euclidean ray, x — n - By this, fl31|) is extended onto 
the whole cut disk. The resulting ratio is equal to J l,m (n, p)/(p n /{' m ); some numbers are 
given in Table 3 to illustrate the effect. While Table 2 shows what improvement of (|30| ) 
is achieved thanks to ([31]) and the Theorem 1, Table 3 shows how strong the standard 
assumption of an overall validity of (|3TD is. Out of the three bounds, Mp n l l { m is the most 
restrictive one, while our result MJ l ' m (n, p) based on the Theorem 1 is looser (as is shown 
in Table 3), and M/{' m based on ( j30|) is the loosest (as is seen from Table 2). Thus, 

M p n l[' m < M J l ' m {n, p)<M l[' m . (37) 
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The efficiency of our bound depends on the value of p, the ratio of 1/d (the lowest value 
of \s\ at which the bound ( |5DD is supposed to hold) to (the radius of the integration 



contour of ( jTT|) in Fig. 1). Clearly, the problem has sense only for p < 1, when the 
integration contour lies in the analyticity region. If the input bounds ( PU| ) and (31) are 
valid down to very low energies (i.e., if 1/d is small with respect to m^), p will be small 
and the right-hand side of (|32|) will become small as well, thanks to the factor p nx ^ . This 
improving factor is, however, strongly angle-dependent in the complex s plane and will 
eventually, near the cut, rise to unity, thereby raising the bound fl32|) back to the starting 
inequality ( PU|) at the points of the cut. 

These results will be modified when the logarithmic dependence of the coefficient 
functions ak(q) is taken into account. 



5 Concluding remarks 

It has been our aim to propose a framework allowing a quantitative estimate of the trun- 
cation error of the operator-product expansion away from euclidean region. In considering 
the evaluation of the contour integrals of the type (|9], [U]) or (|T2|), we have proposed two 
different sets of model assumptions to estimate the influence of the cut s > on the 
truncation error along different rays in the complex s plane. In either case, the starting 
relation is the inequality @ for negative s, s < — b < 0. When combined with analyticity, 
(|) can be extended into the complex s plane, but additional assumptions are necessary. 
We have cosidered two sets of them, (|13"D, (|14D (see Case I of section 4), and (|30|) (see Case 
II). The assumptions are rather strong in both cases, but are weaker than a straightfor- 
ward extension of the inequality (||) into the s plane. Also the resulting estimates (([16]), 



(|l9l) and, respectively, (p2|) ) are considerably looser than such a straightforward extension. 
Some examples to illustrate this are in Table 3. 

This result indicates, in view of the fact that our bounds can be saturated, that in 
applying an operator product expansion one should not mechanically extend ([]) into the 
complex s plane without special physical justification. 

In either case, the resulting bounds exhibit a pronounced angle dependence in the 
s plane, becoming worse and worse with increasing deflection from the euclidean region 
down to the cut. As a consequence, the integrated bounds on (|| [□]) or (|i~2| ) are looser than 
those based on the a priori assumption of angle independence. This result can be under- 
stood as a warning that conventional truncation error estimates used in phenomenology 
are perhaps too optimistic. 

As our results are valid for certain classes of functions, a reduction of the class con- 
sidered to a smaller one could yield a more restrictive upper bound. It is a challenge for 
physics to look for physically motivated class reductions. 

Our results do not represent the complete solution to the problem. As mentioned 
in the Introduction, it was our aim to propose a model scheme that (i) would have 
precise mathematical meaning, (ii) would be free from the a priori assumption that the 
discontinuity along the cut has no influence on the truncation error along a general ray 
and, simultaneously, (in) would keep the model possibly close to real situations. Whereas 
our model scheme satisfactorily meets the conditions (i) and (ii), it does not do full justice 
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to the requirement ( Hi) . We have made a step out of the schemes based on the a priori 
assumption that the truncation error is the same along all rays \s\ — > oo passing through 
the origin. Our approach has been based on a plausible, but still crude picture of the 
operator product expansion. Further refinement is necessary. Next step should include 
allowance for the logarithmic g-dependence of the coefficients a>k{q), and also a discussion 
of current QCD models ||14|| , [|l8|j . An improvement of the integrated error estimate 
may be reached by introducing a subtracted dispersion relation and reversing the order 
of the s-integration and the t- integration (inclusive processes). Work along these lines is 
in progress. 
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Smirnov and J. Stern for stimulating discussions, and M. Flato for drawing our attention 
to Ref. H]. One of us (J.F.) is indebted to A. de Riijula and the CERN Theory Division 
for hospitality. We acknowledge the support of the grants Nos. |gaav-al010711| , MSMT- 
VS96086, and GACR-202/96/1616. 

A Appendix: Theorem A and its proof 

Let C be the plane of complex numbers and C(d, a) = {z G C : z = re lip , \cp\ < a, < r < 
d} where d > and < a < ir are given (see Fig 2). 

Theorem 2. Let f(z) be holomorphic in C(d, a) and fulfilling 

\f{z)\<M (38) 

for z G C(d, a) and 

\f(x)\<Nx n (39) 
for < x < d, where M and iV are positive constants. Denote d = ^ M/N. Then 

|/(*) | < \z\ n/2 VNM (40) 

for z G C(d,a/2), \z\ < mm(do,d) and 

\f{z)\<M for z eC(d,a/2), mm(d , d) < \z\ < d. (41) 

Remark 1. Repeating the statement of the Theorem we conclude that under the condi- 
tions of the Theorem we have: for every nonnegative integer s 

\f(z)\ <N s \z\ n/23 (42) 

for z G C(d, a(l — 2 _s )), \z\ < min(cfo, d) and 

1/(2)1 <M for z EC(d,a(l-2- s )), mm(d , d) < \z\ < d, (43) 

where 

N s = M(N/M) 1/2S (44) 

is valid. This means that the remainder R n {z) (^]) is bounded in the whole region C(d, a) 
but the estimates are bad near the boundary. 
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Remark 2. By a reasoning similar to that used in the Proof of the Theorem the 
following statement can be proved. If a function f(z) fulfils the assumptions of the 
Theorem, then 

\f(z)\<N s \z\^- 2 - s ) (45) 
for z G C(d,a/2 S ), \z\ < mm(d ,d) and 

\f{z)\<M for z eC(d ,a/2 s ), mm(d , d) < \z\ < d. (46) 

is valid where 

N s = M(N/M) l ~ 1/2a . (47) 

This means that, in a small angle, the estimate can be improved. 
Combining these two remarks we have 

Corollary. If a function f(z) fulfils the assumptions of the Theorem, then for every 
£, < £ < 1, the following inequalities 

|/(2)| < N^\z\ n ^ for zeC(d,ag), \ z \ < min(d , d) 

\f(z)\<M for zeC(d,a£), mm(d , d) < \z\ < d (48) 

are valid where 

iV 5 = M(N/M) l -t. (49) 
Proof of the Theorem. Denote 

h{z) = In \f(z)\ -InM. (50) 

The symbol fi(z) is understood as the harmonic function fi(x, y), where z = x + iy. The 
symbol fi(x) is understood as fi(x, 0). f\ is a harmonic function in the region C(d, a) —Af 
where Af — {z : f(z) = 0} is a countable set without condensation points in C(d, a). The 
function /i(.) fulfils 

f\(x) < min(nlnx + lnA^ - lnM,0) (51) 
for real x in (0, d) — Af. Further, 

h{z) < (52) 
in C(d, a) with the exception of countably many isolated points. Define 

G(x) = min(nlna; + lnAr-lnM,0). (53) 

The estimate G(x) of fi(x) is nonpositive for < x < d. 

Now we introduce the notation V = {z e C : z = re llp , < ip < a,0 < r < d}, and 
define the function gi(z) which is 

(a) harmonic, nonpositive and maximal in the region V(d, a) and continuously extensible 
on {x : < x < d}, a part of the boundary ofV(d, a), with exception of countably many 
isolated points, and 

(b) fulfilling 

gi (x) < G(x) for 0<x<d. (54) 
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Such function really exists (see Appendix B). It follows that the function fulfils 

gi(x) = G(x) for < x < d, (55) 
gi(xexp(ia)) = for < x < d, 
gi(d exp(i(p)) = for < <p < a. 

In a similar way, we define a symmetric harmonic function #2(2) fulfilling (a) and the 
inequality 

g 2 (xe ia ) < G(x) for 0<x<d (56) 
with exception of countably many isolated points. This function fulfils 

g 2 (xe ia ) = G{x) for 0<x<d (57) 
#2 (x) = for < x < d 
g 2 (de 1(p ) = for 0<p<a. 

The sum g(z) = g\{z) + g 2 (z) certainly fulfils the inequality (the gi are maximal) 

g(z) < G(\z\) (58) 

so that the symmetry of the functions gi yields 

9l (xe ia / 2 )=g 2 (xe ia / 2 )<G(x)/2. (59) 

If we compare the functions fi and g\ we obtain (see Appendix B) 

fi{z)<gi{z) in V{d,a)-M. (60) 

We have 

gi(xe ia / 2 ) < G(x)/2. (61) 

Since G(x) is nonpositive we have 

9i(x) < G(x)/2. (62) 

These conditions together with (|60| ) and the maximality yields 

fx(z)<G(\z\)/2 in V(d,a/2) (63) 

under the condition that we define fi(z) = —00 for z G M . Using the definitions of the 
functions f\ we obtain the statement of the Theorem. 

B Appendix: Proof of two statements 

(i) Proof of maximality. The region T>(d, a) can be conformally mapped on the unit 
disk. Denote T and 7(.) the image of the interval (0, d] and the image of the function 
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G(x), respectively. Due to the theorem |L9| (Chap. 3, page 33-34, item (v))[]the maximal 
harmonic function fulfilling the condition (a) is given by the Poisson formula 



^L^ \-2\ Z \col(J g ^-^ + \z\^ - (64) 

(ii) Proof of (|60|) . Let {d n }, {a n } and {p n } be sequences fulfilling the inequalities 
d n < d, a n < a and p n > and converging to d, a and 0, respectively. Denote V n ^ m = 
V(p m ,d n ,a n ) = V(d n ,a n ) — V(p m ,a n ). Assume that the sequences are chosen so that 
(YDn,m H Af C (p m ,d n ], where dV n/m denotes the boundary of T> n<m . This means that 
P n>m contains only a finite number of points from M . T> n m can be conformally mapped 
on the unit disc (denoted by U). Further denote f n,m and h n,m the image of f\ and f/M 
respectively. Then certainly f n ' m = In |/i n,m |. The function h n ' m has only a finite number 
of zero points in U. Let {zi, z 2 , z s } be the zero points with multiplicities {ni, tii, n s }. 
Since h nm /~Q.{z — Zi) ni is a nonzero holomorphic function in U the function 



F n > m (z) = f n > m (z)-Y, n M^—^\ (65) 
i=i 1 ZZi 



is a harmonic function in U. Since the functions In \ (z — Zi)/(1 — zz~i)\ vanish on dU and 
f n,m (z) are nonpositive, the functions F n ' m (z) are nonpositive, too. Due to the theorem 
from |19! there exists a nonpositive measure /x on &D such that 



(66) 

Let r n ' m and 7 n,m be the image of the set (p m , d n ) and of the function G(.), respectively. 
Using again the fact that the functions \n\(z — Zj)/(1 — zzj)\ are zero on dU we have 
F n ' m (z) = f n ' m (z) < 7 n ' m (z) on T n ' m and 

F n > m (z) < — [ 7 n ' m (^) rn — , N — ^rd^. (67) 

Z7T Jr™' m 1 — 2\z\ cos(arg2 — ip) + \z\ 2 

Since the functions In \(z—Zi)/(l—zzi)\, k = 1, 2, s, are nonpositive in U (with exception 
of the point z^), we obtain 

f n ' m (w) < q n > m {w) (68) 

in U, where the q n,m (w) are defined by the last integral. Denote f3 n,m the preimage of the 
function q n,m (w). Certainly we have 

fi(z)<p n ' m (z) for zeV (69) 

Due to (i) the functions (3 n,m {.) fulfil the condition (a) of the Appendix A in the region 
V n ^ m and the condition f3 n ' m (x) < G(x) for x G (p m ,d n }. Further, the functions f3 n ' m (x) 

4 A harmonic function / in the open unit disc is the Poisson integral of a finite positive Baire measure 
if and only if / is non-negative. 
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fulfil 



(3^( x ) = G{x) for p m < x < d n , (70) 
/3 n ' m (xexp[ia n }) = for p m < x < d n , 
f3 n ' m (d n exp[iip}) = for < ip < a n , 
/r< m (p m expM) = for 0< V <a n . 

Since the functions /3 n,m (.) are maximal they form a nonincreasing sequence for m — > oo 
and fixed n in V(d n ,a n ), and are bounded from below by the maximal function (3 n in 
V(d n ,a n ), which fulfils (a) in V(d n ,a n ) and (3 n (x) < G(x) for x G (0, d n ]. Denote 
u n (z) = inf m f3 n ' m (z) with fixed n; z/™(.) is a harmonic function in T>(d n ,a n ), and since 
f3 n (x) is maximal we have z/ n (z) = f3 n (z). Certainly we have 

/i(z)</T(z) (71) 

for all n. This procedure can be repeated for n —>■ oo and the inequality fl60|) is proved. 



C Appendix: A function saturating the bound ( |40| ) 



The following example shows that the exponent in ([40]) cannot be improved. 
Example. Let /(.) be the function 

/(z)=exp{ia[ln(v^e-/ 2 )] 2 } (72) 

where a is a positive constant. The function is holomorphic in C(oo,7r) and is bounded 
in C(l,7r). The function f(z) can be rewritten 



f( z ) = \yfl |2a(V2-ar gv ^) eX p{ia [(In ^1 ) 2 - (yr/2 - arg ^ ) 2 ] } (73) 

such that 

|/(2)| = | z |a(ir-arg*)/2_ ( ?4 ) 

This yields that /(x) = x a7r / 2 along the positive real axis and converges to zero as 
| z |a(7r-o)/2 a ] on g ra y ^ = . arg^ = a ; Rez > 0}. Compare this with Corollary. 
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